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Abstract 

The energy spectrum of neutrino-induced upward-going muons in MACRO was analysed in terms of special 
relativity principles violating effects, keeping standard mass-induced atmospheric neutrino oscillations as the 

dominant source of v M <-> v T transitions. The data disfavour these exotic possibilities even at a sub-dominant 

level, and stringent 90% C.L. limits are placed on the Lorentz invariance violation parameters. These limits can also 
be re-interpreted as upper bounds on the parameters describing violation of the Equivalence Principle. 


1 Introduction 

The phenomenon of neutrino flavour oscillations, induced by flavour-mass eigenstate mixing, is now 
considered a solid explanation of both solar and atmospheric neutrino data [1, 2, 3, 4, 5]. In particular, 
the deficit of atmospheric muon neutrinos is well described by the two flavour *—> v T oscillations, 

with maximal mixing between flavour and mass eigenstates [3, 4]. This solution is strongly favoured over 
alternative solutions, like < —> v sterile oscillations [6, 7], <—> v e oscillations [4, 5] and other exotic 

possibilities, such as violation of relativity principles [8, 9, 10, 11, 12], neutrino decoherence [8, 13] and 
neutrino radiative decays [13, 14]. 

However, these alternative mechanisms have been considered under the hypothesis that each one solely 
accounts for the observed effects. It means that they are considered independent and mutually exclusive. 
Here we analyse the possibility of a mixed scenario: one mechanism, the mass-induced flavour oscillations, 
is considered dominant and a second mechanism is included in competition with the former. We studied, 
as sub-dominant mechanism, neutrino flavour transitions induced by violations of relativity principles, i.e. 
violation of the Lorentz invariance (VLI) or of the equivalence principle (VEP). 

In this mixed scenario, we assume that neutrinos can be described in terms of three distinct bases: 
flavour eigenstates, mass eigenstates and velocity eigenstates, the latter being characterised by different 
maximum attainable velocities (MAVs), and consider that only two families contribute to the atmospheric 
neutrino oscillations. Thus, we may write 

K) = K) cos (9* + \ v z) sin0i m 

\v T ) = -\vQ sinOi + |i/$) cosOi 

In Eq. (1) the index * = m or i = v for mass and velocity eigenstates, respectively. When both mass- 
induced and VLI-induced neutrino oscillations are considered simultaneously, the v,, survival probability 
can be expressed as [10, 11, 12] 

P^^ = 1 - sin 2 20 sin 2 Q (2) 

where the global mixing angle 0 and the term Q are given by: 

20 = arctan(ai/a 2 ) . . 

fl = sj (a? + a|) 
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Figure 1: Energy dependence of the <—> v T oscillation probability for mass induced oscillations alone 
(continuous curves), and mass-induced + VLI oscillations for Av = 2 • 10 -25 and different values of the 
sin 2 20 v parameter (dashed curves). The neutrino pathlength was fixed at L = 10 4 km. 


The terms a± and 02 in Eq. (3) contain the relevant physical information 

ai = 1.27|A?n 2 sm20 m L/E + 2 ■ 10 18 Awsin 20 v LEe iri \ 

ai = 1.27 (Am 2 cos 20mL/E + 2 • 10 18 Aw cos 26 V LE') , ' 4 ' 

where the muon neutrino pathlength L is expressed in km, the neutrino energy E in GeV and the oscillation 
parameters Am 2 = m 2 m — m 2 m and Aw = w„« — w„j are in eV 2 and c units, respectively. The unconstrained 
phase i] refers to the connection between the mass and velocity eigenstates. 

Note that in the limit Av = 0 (Am 2 = 0) we recover pure mass-induced (VLI-induced) oscillations and 
the functional form of the oscillation probabilities exhibits an L/E v (L ■ E v ) dependence on the neutrino 
energies and pathlengths. Moreover, in the limiting cases of pure oscillations, the survival probabilities 
do not depend on the sign of the mixing angle and/or on the sign of the Av and Am 2 parameters; this 
is not so in the case of mixed oscillations, as can be seen from Eq. (4), where the relative sign between 
the mass-induced and VLI-induced oscillation terms is important. The whole domain of variability of the 
parameters can be accessed with the requirements A?n 2 > 0, 0 < 0 m < 7 r/ 2 , Aw > 0 and — 7 r /4 <6 v < 7 t/ 4 . 

The same formalism also applies to violation of the equivalence principle, after substituting Aw/2 with 
the adimensional product |</>| A'y; A 7 is the difference of the coupling constants for neutrinos of different 
types to the gravitational potential </ [15]. 

As shown in [11, 12, 16], the most sensitive tests of VLI can be made by analysing the high energy tail 
of atmospheric neutrinos at large pathlength values. As an example, Fig. 1 shows the energy dependence 
of the Vp <—► v T oscillation probability as a function of the neutrino energy, for neutrino mass-induced 
oscillations alone and for both mass and VLI-induced oscillations for Av = 2 ■ 10 ~ 25 and different values 
of the sin 2 2 6 V parameter. Note the large sensitivity for large neutrino energies and large mixing angles. 
Given the very small neutrino mass $ 1 eV), neutrinos with energies larger than 100 GeV are extremely 
relativistic, with Lorentz 7 factors larger than 10 11 . 


2 























7 


• MACRO data 
— Bartol flux 96 
Honda flux 2001 
Bartol flux 96 oscillated 



cos 0 


Figure 2: Comparison of the MACRO upward-throughgoing muons (black circles) with the predictions 
of the Bartol96 and of the Honda2001 MC oscillated and non oscillated fluxes (oscillation parameters 
A to 2 = 2.3 ■ 10" 3 eV 2 and sin 2 2 9 m = 1). 


In the next section we summarise the main results of the global analysis of MACRO neutrino data [3], 
which used the information from i) the zenith angle distribution of upward-throughgoing muons, ii) from 
the energy distribution of a sub-sample of them and iii) from the low energy neutrino data. 

In Sect. 3, we consider the sub-sample of events in point ii) and re-analyse their energy distribution in 
terms of mass-induced + VLI-induced mixed oscillations [17]. Since we consider mass-induced oscillations 
as the dominant source of neutrino flavour transitions, we fix the parameter Am 2 and sin 2 2 9 m at the values 
obtained from the global analysis and leave the VLI parameters free to vary. Finally, the constrains on Am 2 
and sin 2 2 9 m parameters are partially relaxed and a multidimensional analysis is performed. 


2 MACRO results on mass-induced atmospheric neutrino 
oscillations 

MACRO [18] was a multipurpose large area detector (~ 10000 m 2 sr acceptance for an isotropic flux) located 
in the Gran Sasso underground Lab, shielded by a minimum rock overburden of 3150 hg/cm 2 . The detector 
had global dimensions of 76.6 x 12 x 9.3 m 3 and used limited streamer tubes and scintillation counters to 
detect muons, u^s were detected via charged current interactions + N —> fi + X ; upward-throughgoing 
muons were identified with the streamer tube system (for tracking) and the scintillator system (for tirne- 
of-flight measurement). Early results on atmospheric neutrinos were published in [19] and in [2] for the 
upward-throughgoing muon sample and in [20] for the low energy semi-contained and upgoing-stopping 
muon events. Matter effects in the *—> vsterile channel were presented in [6] and a global analysis of all 
neutrino data in [3]. 

In the following we shall briefly recall the analyses performed on different samples of upgoing muon data, 
whose combination yielded the mass-induced atmospheric neutrino oscillation parameters: Am 2 = 2.3 • 10~ 3 
eV 2 and sin 2 2 9 m = 1. 

The first method, (already used in early analyses) is based on the study of the zenith angle distribution 
of upward-throughgoing muons. Fig. 2 shows MACRO data (black circles) compared with two Monte 
Carlo (MC) predictions: the Bartol96 [21] flux with and without oscillations (the dashed and solid lines 
respectively) and the Honda2001 flux [22], the lower curves; the new FLUKA MC predictions agree perfectly 
with the new Honda [25]. The shapes of the angular distributions are very similar, but the absolute flux 
values predicted by the two simulations differ by about 25%; more specifically, if we consider the flux 
extrapolation towards the horizontal direction where no oscillations are expected, we find that our data are 
in good agreement with the Bartol96 MC. 
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Figure 3: (left) Ratio Data/MC n o osc as a function of the estimated L/E v for the upthroughgoing muon 
sample (black points). The solid line is the MC expectation assuming Am 2 = 2.3-10 -3 eV 2 and sin 2 2 8 = 1. 
The last point (empty circle) is obtained from the IU sample, (right) Interpolated 90% C.L. contour plots 
of the allowed regions in the sin 2 29 — A77? 2 plane for the MACRO data using only the ratios R\,R 2 , R 3 
(continuous line) and adding also the information on the absolute values R 4 , R$ (dotted line). 


For a subsample of 300 upward-throughgoing muons MACRO estimated the muon energies via multiple 
Coulomb scattering in the 7 horizontal rock absorbers in the lower apparatus [23, 24]. The energy estimate 
was obtained using the streamer tubes in drift mode, which allowed to considerably improve the spatial 
resolution of the detector (~ 3 mm). The overall neutrino energy resolution was of the order of 100%, mainly 
dominated by muon energy losses in the rock below the detector (note that (E^) ~ 0.4 {E v )). Upgoing muon 
neutrinos of this sample have large zenith angles (> 120 °) and the median value of neutrino pathlengths is 
slightly larger than 10000 km. 

Fig. 3a shows the ratio data/MC as a function of the estimated L/E v for the upthroughgoing muon 
sample. The black circles are data/Bartol96 MC (assuming no oscillations); the solid line is the oscillated 
MC prediction for Am 2 = 2.3 • 10 -3 eV 2 and sin 2 2 9 m = 1. The shaded region represents the simulation 
uncertainties. The last point (empty circle) is obtained from semicontained upward going muons. 

In order to reduce the effects of systematic uncertainties in the MC absolute fluxes we used the following 
three independent ratios [3]: 

(i) High Energy data: vertical/horizontal ratio, R\ = N vert /Nhor 

(ii) High Energy data: low energy/high energy ratio, R 2 = Ni ow /N high 

(iii) Low Energy data: R3 = (. Data/MC)w/ {Data/MC) id+ugs 

With these ratios, the no oscillation hypothesis has a probability P ~ 3 • 10 ~ 7 and is ruled out at the ~ 5cr 
level. By fitting the 3 ratios to the <—» v T oscillation formulae we obtain sin 2 2 9 m = 1, Am 2 = 2.3 • 10 -3 
eV 2 and the allowed region indicated by the solid line in Fig. 3b. 

Using the Bartol96 flux we may add the information on the absolute fluxes of the 

(iv) High Energy data (systematic scale error of ~ 17%): JS 4 = N meas /NBartoi96- 

(v) Low Energy semicontained muons, with a systematic scale error of 21%, R 5 ~ N meas /NBartoi96- 

These informations reduce the area of the allowed region (dashed line in Fig. 3b), do not change the best 
fit values and bring the oscillation significance to ~ 6 a. 
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3 Exotic contributions to mass-induced oscillations 


In order to analyse the MACRO data in terms of VLI, we used the aforementioned subsample of 300 
upward-throughgoing muons whose energies were estimated via Multiple Coulomb Scattering. 

We used two independent and complementary analyses: one based on the y 2 criterion and the Feldman 
and Cousins prescription [26], and a second one based on the maximum likelihood technique [17]. 


3.1 x 2 Analysis 

Following the analysis in Ref. [24], we selected a low and a high energy sample by requiring that the 
reconstructed neutrino energy E^ ec should be E™ 0 < 30 GeV and E^ ec > 130 GeV. The number of events 
surviving these cuts is Ni ow = 49 and Nhigh = 58, respectively; their median energies, estimated via Monte 
Carlo, are 13 GeV and 204 GeV (assuming mass-induced oscillations). 

The analysis then proceeds by fixing the neutrino mass oscillation parameters at the values obtained 
with the global analysis of all MACRO neutrino data [3]: Am 2 = 2.3 • 10 -3 eV 2 , sin 2 20 m = 1. The factor e zr> 
for the moment is assumed to be real (77 = 0 or ir). Then, we scanned the plane of the two free parameters 
(Aw, 9 V ) using the function 


X 


2 



Nj - aN^ c (Av, 9 V ; Am 2 ,9 m ) \ 

J 


(5) 


where N^ IC is the number of events predicted by Monte Carlo, a is a constant which normalises the 
number of Monte Carlo events to the number of observed events and a is the overall error comprehensive 
of statistical and systematic uncertainties [17]. 

We used the Monte Carlo simulation described in [24] with different neutrino fluxes in input 
[21, 22, 27, 28]. The largest relative percentage difference of the extreme values of the MC expected ratio 
Niow/Nhigh is 13%. However, in the evaluation of the systematic error, the main sources of uncertainties for 
this ratio (namely the primary cosmic ray spectral index and neutrino cross sections) have been separately 
estimated and their effects added in quadrature (see [24] for details): in this work, we use a conservative 
16% theoretical systematic error on the ratio Ni ow /Nhigh- The experimental systematic error on the ratio 
is estimated to be 6 %. In the following, we show the results obtained with the computation in [28]. 

The inclusion of the VLI effect does not improve the y 2 in any point of the (Av, 9 V ) plane, compared 
to mass-induced oscillations stand-alone, and proper upper limits on VLI parameters were obtained. The 
90% C.L. limits on Aw and 0 V , computed with the Feldman and Cousins prescription [26], are shown by the 
dashed line in Fig. 4. 

The energy cuts described above (the same used in Ref. [24]), were optimised for mass-induced neutrino 
oscillations. In order to maximise the sensitivity of the analysis for VLI induced oscillations, we performed 
a blind analysis, based only on Monte Carlo events, to determine the energy cuts which yield the best 
performances. The results of this study suggest the cuts E™ c < 28 GeV and El ec > 142 GeV; with these 
cuts the number of events in the real data are N' low = 44 events and N' high = 35 events. The limits obtained 
with this selection are shown in Fig. 4 by the continuous line. As expected, the limits are now more stringent 
than for the previous choice. 

In order to understand the dependence of this result with respect to the choice of the Am 2 parameter, 
we varied the Am 2 values around the best-fit point. We found that a variation of Am 2 of ±30% moves 
up/down the upper limit of VLI parameters by at most a factor 2. 

Finally, we computed the limit on Aw marginalized with respect to all the other parameters left free to 
change inside the intervals: Am 2 = Am 2 ± 30%, 9 m = 9 m ± 20%, — 7 r /4 < 9 V < 7 r/4 and any value of the 
phase 77 . We obtained the 90% C.L. upper limit |Aw| < 3 • 10~ 25 . 


3.2 Maximum Likelihood Analysis 

A different and complementary analysis of VLI contributions to the atmospheric neutrino oscillations was 
made on the MACRO muon data corresponding to parent neutrino energies in the range 25 GeV < E v < 75 
GeV. This energy region is characterised by the best energy reconstruction, and the number of muons 
satisfying this selection is 106. These events are outside the energy ranges used in the analysis discussed 
in Sect. 3.1, thus the expected sensitivity to VLI (or VEP) contributions to the atmospheric neutrino 
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Figure 4: 90% C.L. upper limits on the Lorentz invariance violation parameter Av versus sin 2 9 V . Standard 
mass-induced oscillations are assumed in the two-flavour *—> v T approximation, with Am 2 = 2.3 • 10 -3 
eV 2 and 9 m = 7 t/4 . The dashed line shows the limit obtained with the same selection criteria of Ref. [24] 
to define the low and high energy samples; the continuous line is the final result obtained with the selection 
criteria optimised for the present analysis (see text). 


oscillations should be lower; on the other hand, the maximum likelihood technique (MLT) has the advantage 
to exploit the information event-by-event (it is a bin-free approach). 

Given a specific hypothesis, MLT allows to determine the set of parameters a of the problem (in our 
case a = (Am 2 ,9 m ,Av,9 v )) that maximises the probability of the realization of the actual measurements 
x (here x = (E,L)); this can be accomplished by minimising the negative log-likelihood function: 

n 

£=- 2 ^ 1 n/(x i; a) ( 6 ) 

2=1 

where the sum is over the number of observed events, and /(xp a) is (at least proportional to) the probability 
of realization of a given event, which in our case is: 

/(x; a) =/C(a) • V M C (x) -P^-^(x;a) (7) 

In Eq. (7), /C(a) is a normalisation factor meant to ensure that the integral of /(x; a) over the observables x 
space is independent on the parameters a (otherwise C would not converge), 7Vfc( x ) is the MC probability 
to observe the event x in the no-oscillation hypothesis and P vv (x; a) is the survival probability given 
by Eq. (2). 

MLT has a drawback: C is not a true goodness of fit estimator, as the more popular \ 2 we used in the 
analysis in Sect. 3.1; at the same time, it has the advantage to be very effective when the experimental 
statistics is limited 1 . 

We have chosen different fixed values of the Am 2 and sin 2 2 9 m mass-oscillation parameters on the 90% 
C.L. border reported in Ref. [3] (the solid curve in Fig. 3b) and found the relative Av and sin 2 2 9 V that 

1 We tested the MLT on the 106 events sub-sample of the data in the hypothesis of only mass-induced neutrino oscillations, 

obtaining results compatible with those reported in [3], based on the full MACRO statistics. 
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Figure 5: 90% C.L. upper and lower limits on the Av/2 parameter, versus the Am 2 parameter varying 
inside the 90% C.L. MACRO global result [3]. 


minimise Eq. (6). Fig. 5 shows the 90% C.L. limits of the VLI parameter Av/2 versus the assumed Am 2 
values. 

The limit shown in Fig. 5 was obtained as a convolution of the 90% C.L. upper limits of |Ai>/2| 
corresponding to each chosen Am 2 value. 

4 Conclusions 

We have searched for “exotic” contributions to standard mass-induced atmospheric neutrino oscillations 
arising from a possible violation of Lorentz invariance. We used a subsample of MACRO muon events for 
which an energy measurement was made via multiple Coulomb scattering. The inclusion of VLI effects does 
not improve the fit to the data, and we conclude that these effects are disfavoured even at the sub-dominant 
level [17]. 

Two different and complementary analyses were performed on the data, both of them yielding compatible 
upper limits for the VLI contribution. 

The first approach uses two sub-sets of events referred to as the low energy and the high energy samples. 
The mass neutrino oscillation parameters have been fixed to the values determined in Ref. [3], and we 
mapped the evolution of the y 2 estimator in the plane of the VLI parameters, Av and sin 2 26% No y 2 
significant improvement was found, so we applied the Feldman Cousins method to determine 90% C.L. 
upper limits on the VLI parameters: |Ac| < 6 • 10~ 24 at sin2(9„ = 0 and |Ai>| < 2.5 4- 5 • 10~ 26 at 
sin2d„ = ±1, see Fig. 4. In terms of the parameter Av alone (marginalization with respect to all the other 
parameters), the VLI parameter bound is (at 90% C.L.) |Au| < 3 • 10 -25 . 

These results may be reinterpreted in terms of 90% C.L. limits of parameters connected with violation 
of the equivalence principle, giving the limit |</>A , y| < 1.5 • 10~ 25 . 

The second approach exploits the information contained in a data sub-set characterised by intermediate 
muon energies. It is based on the maximum likelihood technique, and considers the mass neutrino oscillation 
parameters inside the 90% border of the global result [3]. The obtained 90% C.L. upper and lower limits 
on the Av VLI parameter (shown in Fig. 5 versus the assumed Am 2 values) is also around 10~ 25 . 

The two analyses yielded compatible results. 

The limits reported in this paper are comparable to those estimated using SuperKamiokande and K2K 
data [10]. 
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